Proper Scaling of the Anomalous Hall Effect 



o 
o 



o 

CO 



c 

o 
o 



> 
o 

m 
in 

rn 
O 
On 
O 



X 



Yuan Tian, Li Ye, and Xiaofeng JirQ 
Surface Physics Laboratory and Physics Department, Fudan University, Shanghai 200433, China 

(Dated: April 1, 2009) 

Working with epitaxial films of Fe, we succeeded in independent control of different scattering 
processes in the anomalous Hall effect. The result appropriately accounted for the role of phonons, 
thereby clearly exposing the fundamental flaws of the standard plot of the anomalous Hall resistivity 
versus longitudinal resistivity. A new scaling has been thus established that allows an unambiguous 
identification of the intrinsic Berry curvature mechanism as well as the extrinsic skew scattering and 
side-jump mechanisms of the anomalous Hall effect. 
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Shortly after the discovery of the Hall effect, in 1880 
Edwin Hall further observed in ferromagnetic metals an 
additional large contribution besides the ordinary one, 
which is now called the anomalous Hall effect (AHE) - one 
of the most prominent phenomena existing in magnetic 
materials [1] . While the ordinary Hall effect has been well 
understood as a result of the Lorentz force deflecting the 
charge carriers, the mechanism of the AHE has remained 
controversial despite the long history of research, because 
its rich phenomenology defies the standard classification 
methodology, prompting conflicting reports claiming the 
dominance of various processes [2-12]. Recently it again 
attracts great attention because of its natural connection 
to the spin Hall effect and quantum spin Hall effect [13, 
14]. 

In ferromagnets, the transverse resistivity has two con- 
tributions: one is ordinary and is proportional to the ap- 
plied magnetic field; the other is anomalous and is nor- 
mally proportional to the magnetization [8, 9]. It is often 
written as 



Pxy = r H + r a M = p h + p ah 



(1) 



where ro and r a are coefficients that characterize the 
strength of the ordinary and anomalous Hall resistivity 
Ph and pah , respectively. It has long been believed that 
pah should be a function of the longitudinal resistivity 
p xx with well-defined material material dependent pa- 
rameters, i.e., Pah = f(Pxx) ■ 

However, despite the tremendous amount of experi- 
ments, a proper scaling between p a h and p xx has not yet 
been established. In general, p a h could exhibit four types 
of behavior: (a) bp 2 xx (e.g. Fe) [15-18], (b) ap xx + bp xx 
(e.g. Co) [19, 20], (c) bp a xx (e.g. Ni) with 1< a < 2 [21], 
and (d) ap xx (e.g. ultra-pure Ni at low temperature) [22]. 
A unified picture that can explain all the diversifying ex- 
perimental facts is currently absent. 

Theoretically, Karplus and Luttinger first proposed 
that the spin-orbit interaction together with the inter- 
band mixing resulted in an intrinsic anomalous veloc- 
ity in the direction transverse to the electric field [2], 
which gave pi n t oc p xx . This intrinsic contribution 
to the AHE has been recently confirmed in the lan- 



guage of Berry phase [10-12]. However, Smit suggested 
that the skew scattering at impurities was responsible 
for the AHE, which gave p s k oc p xx [3]. Berger fur- 
ther proposed that another impurity-induced mechanism, 
the side-jump, could also give the p S j oc p xx relation 
[5]. In contrast, recent first principles electronic band 
structure calculations based on the Berry phase inter- 
pretation suggested (although implicitly) that it is the 
Karplus-Luttinger intrinsic contribution rather than any 
impurity-induced extrinsic ones that plays the dominant 
role in the AHE [19, 23-25]. Apparently these predic- 
tions are contradictory to each other, while embarrass- 
ingly a direct experimental identification is still lacking. 
In a recent review by Sinitsyn [9], it is shown that the 
total anomalous Hall conductivity (cr a h) consists of five 
different microscopic contributions; they can be further 
divided into the above three categories in terms of the ex- 
perimental identification by transport measurement, i.e., 
o ah = &int+o'sk + Vsj, where <Ti n t is the intrinsic Karplus- 
Luttinger contribution and is irrelevant of any impurity, 
a s k is the extrinsic skew scattering that depends on the 
impurity density, and a S j the generalized extrinsic side- 
jump (including not only the conventional side-jump but 
also the contributions from intrinsic skew scattering and 
the anomalous distribution) that does not depend on the 
impurity density. It is a great challenge to separate the 
roles of intrinsic and extrinsic contributions in experi- 
ment. 

Rising to these challenges we have developed a new 
experimental strategy that goes beyond the existing 
paradigms. Usually the longitudinal resistivity is taken 
as a quantity characteristic only of the material, so tra- 
ditionally besides temperature, p xx was varied only by 
changing the impurity concentration of the material. 
However, such an approach has the unavoidable defect 
that it would in fact modify not only the extrinsic but 
also the intrinsic contributions in the AHE, often com- 
plicating the interpretation of the experimental results. 
Instead, we are tuning the resistivity p xx by varying the 
film thickness of ultrathin layer of Fe, an idea similar 
to the control of the coefficient of viscosity for a mov- 
ing fluid in a thin tube by varying its diameter; the only 
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FIG. 1: (color online) (a) Experimentally measured represen- 
tative p xy vs H curves for 6.5 nm film, from which p a h can be 
obtained, (b) and (c) p a h and p xx as functions of temperature 
for 6.5 nm film, respectively, (d) p a h and p xx for various film 
thicknesses. The dashed blue and solid black lines are fitting 



results with p a n 
respectively. 



b'p xx and p a h = a p xx o + a" p xx r + bp x 



difference is that the former deals with electrons while 
the latter with molecules. It is this novel approach that 
helps appropriately account for the role of phonons in the 
impurity-originated scatterings, resulting in the proper 
scaling of the AHE: a ah = -{a s kV xx0 + b sj a~ x0 )a 2 cx + a vnt 
or pah = {askPxxO + bsjPlxo) ~ G in tp 2 xx ■ Here p xx0 is the 
residual resistivity caused by defects in crystal, a s k and 
b s j are material dependent constants for the skew scatter- 
ing and side-jump, respectively. This new scaling enables 
us to separate out the intrinsic Karplus-Luttinger con- 
tribution from the various impurity-originated extrinsic 
contributions, and further singles out its dominant role 
for the AHE as a xx — ► . We finally develop a unified 
physical picture that can explain all the previous exper- 
imental results. 

Fe films were grown on undoped GaAs(OOl) at 300 
K by molecular beam epitaxy (MBE), and were capped 
with 4 nm thick MgO to prevent oxidation in air. The 
detailed experimental setup was described elsewhere [26, 



27]. The films were patterned into the form of a stan- 
dard Hall bar along [110] with the magnetic field along 
[001]. The transport measurements were carried out in a 
physical property measurement system (Quantum Design 
PPMS-9T system). The magnetoresistance in Fe films at 
5 Tesla is smaller than 0.5%, and the magnetization in 
Fe films thicker than 1 nm has its bulk value and changes 
little within the temperature range of 5 K-320 K. 

The experimental relation between the anomalous Hall 
resistivity p a h and the longitudinal resistivity p xx is es- 
tablished through the following procedure. Fig. 1(a) 
shows several representative sets of p xy vs H curves mea- 
sured at different temperatures between 5 K and 320 K 
for a 6.5 nm thick Fe film. p a h(T) is then obtained as the 
zero field extrapolation of the high field data as shown 
in the figure, and is displayed in Fig. 1(b). p xx (T) was 
measured simultaneously and is shown in Fig. 1(c). From 
Fig. 1(b) and 1(c), p a h versus p xx curve for the 6.5 nm 
Fe film can be deduced and is displayed in Fig. 1(d), 
together with data for other thicknesses varying between 

I. 5 nm and 93 nm. 

This experimentally established p a h = f(pxx) at dif- 
ferent film thicknesses provides an opportunity to unveil 
the phonon contribution to the AHE - a long-standing 
controversial issue [3, 6, 28]. Before going to the de- 
tailed data analysis, we recall first some very basics 
about the anomalous Hall effect. According to Ohm's 
law, there exists a general relation between the anoma- 
lous Hall resistivity and the anomalous Hall conductiv- 
ity: Pah = -Oah/{vl x +<?lh) Or C7 ah = - Pah / \p\ x + p\ h ) , 

which is material independent and is valid for each in- 
dependent mechanism of the AHE. Since a a h << cr xx , 
Pah << Pxx, and a xx = p xx , this can be further simpli- 
fied to: p ah = -<7 a hpl x or a ah = ~PahV 2 xx - It immedi- 
ately follows that different mechanisms are additive both 
in the anomalous Hall conductivity and resistivity, i.e., 
&ah = a j an d Pah = Pi ' an important fact that is 
quite often overlooked. Specifically, a constant <Ji nt there- 
fore p^t = —PintPxx is expected for any given material [2, 

II, 12]; on the other hand, it was generally adopted with- 
out serious justification that the skew scattering and side- 
jump were expressed as p sk = a sk p xx and p sj = b S jp 2 x 
respectively. Following this line of reasoning, it seems 
quite natural why the formula p a h — a<skPxx + bp 2 x has 
been widely used to investigate the AHE and separate 
various contributions (e.g., [18, 19, 25]). But we are go- 
ing to expose the fundamental flaws of this formula in 
the following, then establish a new scaling for the AHE. 

Following the observation in bulk Fe reported in Ref. 
[15, 16], we used b'p 2 x to fit the data in Fig. 1(d), but 
found the significant deviation when the Fe films were 
thinner than 6.5 nm, as clearly seen by the dashed lines 
in the figure. This means that the skew scattering must 
be considered in ultrathin Fe films although it is negligi- 
ble in the bulk material. We therefore use the routinely 
accepted formula p a h = o, sk p xx + bp xx to fit the data, but 
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FIG. 2: (color online). Red, blue, and green dots represent 
respectively the parameters a , a" and b. Solid curves are 
guide to the eye. 



found a s k is not a constant, showing strong thickness 
dependence (not shown here). However since the impu- 
rity density concentration and the boundary condition 
remain the same for different film thickness, it is known 
that p xx decreases as the film thickness increases, but it 
is puzzling why a sk should be thickness dependent or p xx 
dependent (a sk (p xx ) )■ Unlike the experiments with vari- 
able impurity where a s k(p xx ) might happen, here in this 
new approach the fact that a sk is not a constant strongly 
implies the improper scaling of the formula. We start 
to question why the skew scattering should be included 
as p s k = a s kp xx - According to the Matthiessen's rule 
Pxx = Pxxo + PxxT as seen in Fig. 1(c), p sk = a sk p xx 
implies p sk = a sk p xx0 + a sk p xxT , i.e., the phonons would 
contribute (via p xx t) equally to the skew scattering just 
like the defects in crystal (via p xx o). This is certainly a 
very strong assumption and must be verified. In order to 
achieve this we first assume they contribute unequally to 
the skew scattering, i.e., p sk = a' p xx o + a" PxxT , then try 
to determine whether a' and a" are identical by fitting 
the data in Fig. 1(d) with p ah = a'p xx0 + a" p xxT + bp 2 xx 
using the experimentally measured p a h, Pxxo, PxxT, and 
p xx for different film thicknesses, similar to those ob- 
tained for the 6.5 nm film in Fig. 1(b) and 1(c). The 
fits are shown by the solid lines in Fig. 1(d), and the 
fitting parameters of a' , a" and b are presented in Fig. 2. 
It is evident from Fig. 2 that a' and a" are very different 
and actually a" « 0. However, it seems that the same 
puzzle - a' is thickness dependent - still exists; it turns 
out that the two situations are fundamentally different, 
which becomes clear in the following. 

Similarly, the same concern exists for the side-jump as 
well, i.e., why it should be included as p S j = b S jp 2 x in 
the widely used formula, unlike in the intrinsic Karplus- 
Luttinger case where <ii n t is a real constant so that 
Pint — —&intP xx follows directly from the general rela- 



tion between the anomalous Hall conductivity and resis- 
tivity. If the side-jump p S j — b S jp xxQ + b S jTP xxT (it does 
not affect the following results even if there exists cross 
terms of p xx op XX T) is considered in parallel and together 
with the skew scattering p sk = a sk p xxQ + a skT p XX T, then 
without losing the generality what we should verify is 

Pah = {a-skPxxQ + a sk TPxxT) + {bsjP xx Q + b s jTP xxT ) ~ 

dint Pxx ■ If seems unreasonable from the first glance 
to work out so many parameters from Fig. 1(d) alone; 
however, the situation is dramatically changed when 
we regroup the formula in the form of p a h = {&sk + 
b S jPxxo)pxxO + (a s kT + b S jTPxxr)pxxT — &intPxxi which 
has exactly one-to-one correspondence with the formula 
Pah = a'pxxo + a" PxxT + bp 2 xx used to check the skew scat- 
tering above. Apparently we can reinterpreted the same 
Fig. 2 as the following: a" m is a direct experimental 
justification that phonons actually contribute little to the 
overall extrinsic AHE as compared to that of the defects 
in crystal, either the contributions from the skew scatter- 
ing and side-jump cancels each other, or they are both 
negligible; b is nothing but the negative anomalous Hall 
conductivity b — —(Tint, which is fully developed and be- 
comes almost saturated (b = —(Tint ~ 1-1 x 10 3 il _1 crn _1 ) 
as the film thickness reaches 4 nm and above; the thick- 
ness dependent a! in Fig. 2 simply indicates that the 
extrinsic AHE contains contributions not only from the 
skew scattering (a sk p xx o) but also from the side-jump 
(b S jP xx0 ) so that a thickness dependence a' should be 
expected as a' = (a sk + b S jp xx o). Following this line of 
logic, a new scaling of p ah = f(p xx ,p xx o) rather than the 
traditional p a % — f(p xx ) is proposed for the AHE: 

Pah = (a sk p xx0 + b sj p 2 xxQ ) - <Jintpl x (2 A) 

Oah = -(a>sk&xx0 + b sj CJ xx o) (T lx + a int (2-B) 

Here, <Ji n t, a sk and b S j are all material dependent con- 
stants. 

To confirm indeed a' = (a sk + b S jp xx o), i.e., the extrin- 
sic anomalous Hall resistivity p ext — a sk p xx o + b S jp xx0 
or anomalous Hall conductivity a ext — —(a>skC XX Q + 
b s jcr xx0 )a xx as in Eq. 2A and 2B, we plot in Fig. 3 
—Paho vs a xx o as shown by the red dots, using the ex- 
perimental raw data measured at 5 K for different film 
thicknesses. At this temperature, besides p xx o = & xx0 
we also have a ah w a ah0 and a xx w a xx0 , so that Eq. 
2B becomes: -a a ho = a sk a xx0 + (b SJ - a int ), by which 
the experimental data can be well fitted as seen from 
the black line in the figure, meanwhile the correspond- 
ing constants are extracted as: a sk — —3.7 x 10~ 3 and 
(b B j — (Tint) = 1.8 x 10 3 r2 _1 cTO _1 . Recalling the previous 
result —(Tint ~ 1.1 x 10 3 n~ 1 cm~ 1 from Fig. 2, we further- 
more get the side-jump constant b S j w 0.7 x 10 3 fi -1 c7n -1 . 
It is evident that at low temperatures the side-jump b S j 
is not negligible comparing with the Karplus-Luttinger 
intrinsic term —<Ti n u on the other hand, the change of 
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FIG. 3: (color online). The red dots represent the —a a ho vs 
UxxO relation using the experimentally measured raw data at 5 
K for different film thicknesses. The black curve is the fitting 
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FIG. 4: (color online). — a a h vs a%x f° r F e films with various 
thicknesses. The black lines are linear fitting of corresponding 
data. The red dashed line corresponds to the —o a h value 
obtained from iron whisker. 



^sk^xxO as a function of <t xx q for the samples explored 
in this experiment as seen in Fig. 3 can be as large as 
1.0 x 10 3 f2 _1 cm _1 , almost the same magnitude of the 
—Cint value, thus is also not small at all. In addition, 
the negative sign of a s k = —3.7 x 10 -3 indicates that the 
skew scattering contributes to the AHE in Fe in the oppo- 
site direction as the side-jump and the Karplus-Luttinger 
terms do. Therefore in principle it could exceed them and 
becomes dominant at low temperature for samples with 
larger a XX Q, explaining the striking and long puzzled phe- 
nomenon in which the anomalous Hall resistivity of Fe 
would change sign simply as the temperature is lowered 
as observed earlier [16], which was unable to understand 
with the simple p a h = b'p xx term. 

Instead of the data fitting as done in Fig. 1(d) and Fig. 
2, we are going to single out now the Karplus-Luttinger 
intrinsic contribution from the extrinsic contributions in 
a much more straightforward and transparent way. Fig. 
4 shows the a a h versus a xx (T) plot using the experimen- 
tal raw data, each curve corresponding to a specific film 
thickness with variable temperatures between 5 K and 
290 K. This figure contains some important information 
about AHE, which has been hidden too long. First of all, 
the linear relationship between a a h and a xx confirms in 
an elegant way that phonons do contribute little to the 
skew scattering and side-jump in the AHE, as predicted 
by Eq. 2B, otherwise a linear relationship of a a h versus 
a xx would be expected from the widely used but flawed 
formula p a h = a sk p xx + bp xx . Then, as a xx goes to zero 
in the figure, the anomalous Hall conductivity a a h con- 
verges to essentially the same but nonzero value. This 
converged value is nothing but exactly the long searched 
Karplus-Luttinger intrinsic or the Berry curvature con- 
tribution, i.e., a ah = ((Tint + a sk + cr s j) — > a int when 
&xx — * according to Eq. 2B. To flesh out this criti- 



cal point, we believe that as a xx approaching zero (i.e. 
in the high temperature limit for metals), the extrin- 
sic terms (a s k + <J s j) attributed to the impurity scat- 
tering ought to be washed out by the random and in- 
coherent phonon scattering, therefore shrunk to zero as 
seen in Fig. 4; in the meantime, the intrinsic one of (Ji n t 
originated from the electronic structure of the material 
should be the only robust one that is usually less sen- 
sitive to temperature. In addition, this converged value 
of -a int w (1.1 ± 0.1) x lO^^cm- 1 at a 2 xx = for 
films thicker than 4 nm not only agrees very well to 
the (6 = —<Ti n t) result of Fig. 2, but also to that from 
bulk Fe whisker measured at room temperature [16, 24] 
as marked by the red dashed line, which demonstrates 
unambiguously that it is the Karplus-Luttinger intrin- 
sic rather than any extrinsic mechanisms that plays the 
dominant role for the AHE in bulk Fe at room temper- 
ature and higher. It should also be pointed out that for 
films thinner than 4 nm the —<Jmt value decreases as seen 
in Fig. 2 (not shown in Fig. 4), presumably due to the 
finite-size or quantum- well modification (in ultrathin film 
of Fe) to the bulk electronic band structure. 

With the new proper scaling of the anomalous Hall 
effect, Eq. 2 A and 2B, we can now finally unify all the 
aforementioned diverse experimental results in literature. 
First, it is likely that type (a) as defined in the introduc- 
tion corresponds to situations in which either the mate- 
rial dependent parameters a s k and b S j happen to be very 
small, or the measurements were carried out at temper- 
atures where p XX Q <C p xx . Second, we have reanalyzed 
the data of type (b) using Eq. 2A, and found that the 
new scaling can indeed describe better those experimen- 
tal results. Third, type (d) actually belongs to a special 
case of Eq. 2A where the temperature was fixed very 
low (so p xx ~ Pxxo) and the sample was ultra-pure (so 



5 



Pxxo ^ Pxxo), thus the second and third terms in the 
equation are negligible compared to the first, leading to 
an equation of p a h « a s kPxxQ- Fourth, type (c) corre- 
sponds to a nontrivial case of Eq. 2A, where the intrin- 
sic conductivity ai n t itself is sensitive to temperature in 
the range of interest [23, 24]. However, because of the 
absence of the Matthiessen's rule in semiconductors, it 
is not trivial whether the new scaling established here in 
Eq. 2 applies to the extrinsic AHE in magnetic semi- 
conductors as well, although nothing is against it yet. 
Last but not least, if an experiment is carried out in a 
constant longitudinal current mode for an almost perfect 
ferromagnetic crystal at extremely low temperature, the 
anomalous Hall voltage is expected to be essentially zero 
according to Eq. 2 A, but the anomalous Hall current is 
expected to be finite as seen from Eq. 2B. 

We believe that the new result presented here opens 
wide possibilities to manipulate in a controlled way ei- 
ther intrinsic or extrinsic or both effects to meet certain 
application purposes in future spintronics devices. 
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